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Метод параллельно-иерархического 
распознавания изображений 


В статье рассмотрен метод параллельно-иерархического преобразования на основе популяционного 
кодирования и его применение для задач распознавания образов. Параллельно-иерархическое 
преобразование описывается как системная модель для распознавания образов. Приведены 
теоретические сведения, экспериментальные исследования и программная реализация. 


Введение 


Последние достижения нейробиологических исследований показывают, что 
для кодирования сенсорной информации мозг использует подход, который получил 
название популяционного кодирования. При этом информация представляется целой 
популяцией активных нейронов [1]. Это важное свойство подтверждено экспери- 
ментами Д. Спаркса, которые позволили, в результате анализа того, как мозг обезьяны 
руководит движением ее глаз, сделать вывод, что необходимое движение кодируется 
целой популяцией клеток, каждая из которых представляет несколько отличающихся 
от других движений. Движение, которое в результате осуществляется глазом, отвечает 
среднему результату по всем движениям, закодированным активными клетками. 
Эксперименты М. Янга и С. Ямане с височной корой мозга обезьяны показали, что 
популяционное кодирование используется не только для движения глаз, но и для мимики. 

В соответствии с идеей популяционного кодирования образ, который распознается, 
должен сравниваться с некоторым усредненным отображением эталонного образа. 
Построить модель такого усредненного отображения эталонного образа — среднего 
по всем текущим образам (иначе — по всей обучаемой выборке) — можно на основе 
усреднения параметров параллельно-иерархической (ПИ) сети [2]. В представленной 
статье исследуется гипотетическая модель структуризации информации в коре 
головного мозга на основе ПИ сети, обстоятельно изложенной в работе [2], для задач 
распознавания образов. Формирование на ее основе многоэтапной ПИ сети представляет 
процесс последовательного преобразования коррелированных пространственных 
областей и образования декоррелированных во времени элементов нейронной сети 
при переходе из одного устойчивого состояния в другое. Главной особенностью 
предложенного метода является динамика пространственно коррелированного механизма 
преобразования текущих и образования результирующих элементов нейронной сети. 
Такой механизм позволяет представить обработку изображений в нейронной сети 
как процесс параллельно-последовательного преобразования различных составляющих 
образа с учетом их временных характеристик. Физическое содержание входных 
элементов нейронной сети, которые принимают участие в процессе корреляции- 
декорреляции таких, как, например, амплитуда или частота, фаза или энергия сигналов, 
связность или текстура изображений, определяется типом используемого преобразова- 
ния, выбор которого зависит от класса решаемых задач [3]. 
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Описание метода организации ПИ сети 


Анализ образа заключается в последовательном преобразовании совпадающих и 
выявлении (фильтрации) несовпадающих во времени составляющих образа при переходе 
элементов нейронной сети из текущих энергетических состояний с одними простран- 
ственными координатами в состояния с меньшей энергией и другими пространствен- 
ными координатами. Условием перехода составляющих образа на следующий уровень 
является наличие динамики взаимного совпадения промежуточных результатов обра- 
ботки во времени в одноименных каналах нижнего уровня. Результат анализа образа 
формируется из изолированных в пространственно-временной области составляющих об- 
раза [2]. 

Рассмотрим математическую модель параллельного разложения множества 


и = {а} = 1 п [2], [4], применяемого в каждой ветви ПИ сети. 


а» в (’-а”"), (1) 
Е т 1—0 


|: тр 
где а, = 0, В — размерность данного множества, 4’, К =1, К, — элементы подмножеств, 


состоящих из одинаковых элементов, п, — число элементов в А-м подмножестве (то 


р: е ‹ 
есть кратность числа @,), а’- элемент множества {а я выбранный на /-м шаге, 


= Ка =. я 50, 

Воспользуемся идеей популяционного кодирования [1], построив модель 
какого-либо конечного действия, совершаемого по всем текущим действиям. 

Для ПИ сети усредненными параметрами могут быть количество элементов в 
ветви каждого уровня, которые определяются на основе модели вида (1), и значения 
самого элемента. Очевидно, что на уровне ветвей нейронной сети данное конечное 
действие при реализации популяционного кодирования отвечает усредненным 
параметрам этой сети [5]. В этом случае текущий образ, который распознается, будет 
отображаться текущей ПИ сетью и сравниваться с эталонной ПИ сетью с усреднен- 
ными параметрами: 


аа (2) 


Обозначив среднее значение произвольного элемента первого уровня через 
1 о 3 
А:.;, второго уровня — а;.;, третьего уровня — @;.; и т.д. А-го последнего уровня 
р: 
а.;, а также среднее число элементов первого уровня — М, ‚ второго уровня — 
а:.) 


М , , третьего уровня — №, ит.д. А -го уровня — № ‚ , можно сформировать ПИ 


41.) 41.) 41.) 
сеть с усредненными параметрами. 


Синтезированная таким способом структура ПИ сети с усредненными пара- 
метрами показана на рис. 1. 
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Рисунок 1 — Структура ПИ сети [2] с усредненными параметрами 


Текущий образ или исследуемый образ, который обработан ПИ сетью с теку- 


1 2 3 К 
1,1? а, ›› а а и соответствующим числом элементов в 


щими параметрами а рр Ч) 


ветвях каждого уровня — № и № 2 № 8 ...„М „‚ › сравнивается с эталонным образом, 
(7 1.) 1] 1] 
| 2. -3 р 


который обработан ПИ сетью с усредненными параметрами — а:;, @:.;, @г}..... Ч] 


И средним числом элементов В ветвях соответствующего уровня = 


И 


аг.) а.) а1. а1.} 


Если воспользоваться препарированием элементов ПИ сети для каждого 
| 2 3 к 


уровня [3], то от усредненных параметров @;;, @:.;, а, .».... Аг; можно перейти на 


основе трехуровневого кодирования к их отображению бинаризованными препа- 
ратами (—1, 0, +1). Затем определяются массивы разностей элемента со средним 
значением яркостей элементов разложения изображения (или его фрагмента), в 
котором находится элемент изображения: 


вл=а: @ы. (3) 
Для препарирования начального изображения полученные разности сравниваются с 
порогом д ‚то есть 
. если В,,>д 
а) =) -) если В, <-д. (4) 


0, если |. <0д 
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Причем, порог д вычисляется из условия (5): 
(И. А-Ю. ло) — 
‘М ‘И - Мах , (5) 
1 —1 0 

где М”, М”, М” — количества единичных, минус единичных и нулевых препаратов 
при каждом пороге 1 [6]. Как следует из графика на рис. 2, оптимальная величина порога, 
например, для исходного изображения а01_по1.рсх составляет 44. В ходе эксперимента 
были проведены аналогичные операции по оптимизации выбора порога для всех зашум- 
ленных изображений. Результаты выбора оптимальных порогов представлены на рис. 3. 
Гистограмма на рис. 3 показывает, что величина оптимального порога препарирования 


несколько снижается с ростом степени зашумления изображения, причем этот эффект 
проявляется в большей степени для равномерного шумового распределения [6]. 


Выбор оптимального порога: распределение 
препаратов 
(исход.изображение:а1 по1.рсх) 


[-*] 

о 

|= 

5. 

Ы —5 + преп. 
в. —==- - преп. 
о == 

8 0 преп. 
8 произв. 
РУ 

= 

Е 

9 

ы 


1 И 21 31 41 51 61 71 


величина порога 


Рисунок 2 — Функции распределения препаратов (-1,0,1) и их произведения 


Выбор оптимального порога препарирования для 
различных шумовых распределений 
(исход.изображение: а| по1.рсх) 


величина оптим. 
порога препар-я  „ 


Гаусс 
тип шума 
равномер. 


степень зашумления 


Рисунок 3 — Распределение оптимального порога 
для зашумленных изображений 
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Для произвольного усредненного параметра трехуровневый переход может 


0 | 
отображаться тремя видами препаратов: нулевым —а; ‚› положительным — @,;,, 


отрицательным — а, . В этом случае ПИ сеть с числовыми отсчетами преобразуется 
от —1 
‚> Ч; ‚› Ч; ;. Тогда существенно упрощается 
процедура корреляционного сравнения бинаризованных отсчетов текущей и эталонной 
ПИ сетей. 

Для формирования эталонов изображения необходимо провести обучение в 


пределах выборки. Для чего необходимо делать усреднение по элементам ветви 
1 2 3 К 


каждого уровня, то есть формировать усредненные элементы а;;, ау, @лу,.... Чл» 


в ПИ сеть с бинаризованными отсчетами а 


я р а, р в Выполнив вышеуказанные 
действия, можно сформировать ПИ сеть с эталонными параметрами для текущих 
изображений. После чего можно делать ее корреляционное сравнение с ПИ сетью, 
которая использует текущие параметры. (Под ПИ сетью с текущими параметрами 
| ОЕ К 


Ч бь 9 


переходя затем к бинаризованным препаратам а 


понимается ПИ сеть с текущими значениями ее элементов а 


0 1 —1 
1,1? а а, и текущим числом 


переходом к бинаризованным препаратам а 5» Ч, 


элементов в ветвях каждого уровня Мо | М › : М ь ам .) 
1] ] 1] 1] 
Коэффициент корреляции р, между выборками двух случайных дискретных 
величин Х, У с математическими ожиданиями М,,М, и среднеквадратичными 
отклонениями ух, уу определялся как: 


_ СО) 
Ух’Ух 


(6) 


ХУ 


1 Ш 
где -1<р,,<|и сОх)=-У (мм, (у,-м). 


1 
Очень важным обстоятельством является то, что коэффициент корреляции 
можно вычислять не только для каждых двух уровней, но и для двух ПИ сетей, что 
улучшает достоверность результата распознавания. 


Распознавание динамических образов изображений пятен 
лазерных пучков 


В настоящее время в полиграфии, лазерной обработке материалов, локации, опти- 
ческой связи и других областях техники ощущается необходимость более широкого вне- 
дрения оптико-электронных систем с автоматической корректировкой искажений 
формируемого светового излучения. Причинами искажений могут быть дестабилизиру- 
ющее воздействие механических или климатических факторов, нестабильность харак- 
теристик источника излучения, возмущения в оптическом тракте, разъюстировка опти- 
ческих элементов и т.п. Обеспечение приемлемого качества коррекции требует непре- 
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рывного динамического контроля характеристик светового излучения, например про- 
странственного распределения его интенсивности, в том числе оценки отклонения ука- 
занного распределения от исходного или эталонного распределения. 

Рассмотрим распознавание динамических образов на примере изображений 
пятен лазерных пучков. (Рассмотрены исходные сигналы на примере последователь- 
ностей серий изображений протяженных лазерных трасс [7].) 

Решение этой задачи проводилось на основе туннелирования [8] по усредненным 
коэффициентам корреляции при сравнении эталонного образа с отдельным фрагментом 
трассы. Обучение ПИ сети сводилось к выполнению следующих этапов: 

1) для фрагмента трассы найти точку с максимальной яркостью Т 


тах ; 
2) выбирать необходимое количество М значений яркости Т„_,, \=1,2,..., М, 
начиная с максимальной (крайние точки с яркостями 7Т’_, образуют М краевых линий); 


3) для каждого краевого контура определить предварительные центры отдельно 
для каждой координаты: 


х. 5572. . (у). 2% = (И бьз) у 
м=УУ (у) ©) 


х=0 у=0 


ь- 


где и=|,2,..., М, Л (х, у) — функция яркости изображения; 


4) найти усредненные координаты для всего фрагмента трассы (данными для 
усреднения являются координаты центра каждого краевого контура) [8]; 

5) найти эталонный образ (усредненные значения яркостей элементов разложения, 
из которых состоит фрагмент трассы, для обучающей выборки находятся для каждого 
фрагмента трассы); 

6) туннелирование произвести по усредненным коэффициентам корреляции 
при сравнении эталонного образа с отдельным фрагментом трассы (в качестве 
верхней границы туннеля выбирается единица, а в качестве нижней - усредненный 
коэффициент корреляции); 

7) определить «хорошие» или «плохие» изображения (фрагмент трассы считается 
«хорошим» при попадании среднего коэффициента корреляции в туннель эталонного 
образа, иначе — «плохим»); 

8) произвести оптимизацию весов для «хорошего» изображения (8): 


ыы у). («+ 5х), У =— У» У))-(›-+5,), (8) 


М до, М до, 


где и=1,2,..., М, Л (ху) _ яркость, дух, 6, — смещение центров по осям Хи У 


соответственно. 
Смещение центров по осям Х и У определяется по формулам (9) и (10) 
соответственно: 


ь. =, В} (9) 


>.6-х) (10) 
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где х, у — средние значения центра выборки изображений по осям Х и У соответ- 


ственно, х,, у, — координаты центра каждого 1-го «хорошего» изображения. 


Полученные значения х„, у» и есть окончательный результат измерения 


координат энергетических центров изображений пятен лазерных пучков. 


612 


Алгоритм обработки динамических образов приведен на рис. 4. 
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Фрагмент трассы 
“плохой”. 


спешное 
окончание 
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: 


Получение М кривых корре- 
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уровням ПИ-сети. 


> 


Нахождение общего числа 
уровней для Л! фрагментов и 
среднего коэффициента корре- 
ляции по уровням ПИ сети. 


Коэффициенты 
корреляции усреднены? 


нет да 


у 
Построение графика сред- 
них коэффициентов корре- 
ляции по уровням ПИ сети 


у 


Нахождение среднего 
коэффициента корреля- 
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Фрагмент трассы 
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у 
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весов. 


у 
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Рисунок 4 — Алгоритм обработки динамических образов 
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Результаты экспериментальных исследований обработки 
динамических образов изображений 
пятен лазерных пучков 


На рис. 5 представлены образцы эталонных образов и отдельных фрагментов 
протяженных лазерных трасс № 1 — № 10 - изображений пятен по 1000 в каждой из 
трасс [9], которые использовались в ходе экспериментов, а также кривые корреляцион- 
ного сравнения эталонных образов и изображений пятен разных фрагментов трасс по 
уровням ПИ сети. 


Эталонный образ сформирован из 10% фрагментов трассы № 1: 


[3СВ_001.ВМР...$СВ_100.ВМР]; 


Первые 10 фрагментов трассы №1 для сравнения изображений 


по 0 


Рисунок 5 — Образцы образов и фрагменты трасс 


Для фрагментов трассы № 1 строим ПИ сети и получаем кривые корреляци- 
онного сравнения по уровням ПИ сети (рис. 6). 


Коэффициент корреляции 
, , В 
[== 
| 


Ри: 
13579 12 16 20 24 28 32 36 40 44 48 52 56 60 64 68 72 76 80 84 88 92 96 


Номер уровня. 
Рисунок 6 — Пример корреляционного сравнения по уровням ПИ сети 


Для рис. 6 общее число уровней ПИ сети - 18. Далее построим график средних 
по уровням коэффициентов корреляции (рис. 7). 
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Коэффициент корреляции 


Рисунок 7 — График средних по уровням коэффициентов корреляции 


Получили: средний коэффициент корреляции по уровням ПИ сети - 0,41, 
следовательно, туннель находится в пределах - [0,41...1]. 

Определяем «хорошие» фрагменты трассы № 1 и для них проводим оптимизацию 
весов по формуле (8), тем самым вычисляя координаты энергетических центров 
изображений пятен лазерных пучков (рис. 8). 


Значение координат по осям ХиУ 


1 35791216 20 24 28 32 36 40 44 48 52 56 60 64 68 72 76 80 84 88 92 96 
Фрагмент трассы 


Рисунок 8 — График определения координат энергетических центров изображений 
пятен лазерных пучков до и после обучения ПИ сети 


Выводы 


В статье предложен метод организации ПИ сети, использующий идею популяци- 
онного кодирования в искусственной нейронной сети и приближения ее к естествен- 
ным нейронным сетям. Предлагается текущий образ представлять текущей ПИ 
сетью с текущими параметрами и преобразовывать их на основе обобщенного кон- 
турного препарирования [10] в бинарные препараты с последующим корреляционным 
сравнением с эталонным образом эталонной ПИ сети с усредненными параметрами, 
элементами которого являются бинаризованные препараты. В отличие от известных 
структур искусственных нейронных сетей [11], в которых для сравнения используются 
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ненормируемые (абсолютные) критерии подобия, в рассмотренном методе организа- 
ции ПИ сети используется нормируемый критерий. При этом коэффициент корре- 
ляции предложено вычислять не только для каждых двух уровней, но и для двух ПИ 
сетей, что улучшает точность результата распознавания. Результаты машинного 
моделирования, проведенные для примера обработки изображений человеческих лиц, 
показывают полную адекватность полученных теоретических и экспериментальных 
результатов. 

При использовании параллельно-иерархической модели для задач распознавания 
динамических образов обучение сети показывает, что среди «хороших» фрагментов 
трасс оказались примерно 50% фрагментов трасс, для которых была проведена опти- 
мизация весов. Проведенная оптимизация весов значительно увеличивает точность 
определения координат центров изображений пятен фрагментов лазерных пучков. 

Эксперименты показывают, что для изображений пятен лазерных пучков в ре- 
зультате действия разнообразных дестабилизирующих факторов невозможно точно 
измерять координаты их энергетических центров, однако точность можно существенно 
повысить, используя обучение ПИ сети с последующей их классификацией на «хо- 
рошие» и «плохие». 

Сравнительный анализ показывает, что предложенные методы позволяют измерять 
координаты центра лазерных изображений с точностью определения точки привязки 
не более чем 1,2 элемента разложения, что по точности почти на порядок превышают 
известные, например, на основе определения центра масс с помощью моментных 
признаков. Физическое моделирование метода определения координат изображений 
пятен лазерных пучков с применением сигнального процессора типа ТМ5320С5510 с 
тактовой частотой 200 МГц и производительностью 400 млн оп/с показывает, что 
для обработки одного изображения нужно -10,65мс [9], что является целиком доста- 
точным для реализации телевизионной обработки изображений лазерной трассы. 

Полученные результаты могут найти широкое применение в разных системах, для 
которых необходимо осуществлять автоматическую коррекцию искажений лазерного 
источника светового излучения. Это, например, системы оптической навигации с 
использованием лазеров в условиях со сниженной видимостью в атмосфере, системы 
точного ориентирования, которые дают возможность определить пространственное 
положение движущегося объекта в разных метеорологических условиях, лазерные 
навигационные системы. Это также лазерные системы связи и передачи информации, 
локации и дальнометрования, трассирования и т.п. системы, которые испытывают 
влияние атмосферы на параметры лазерного луча, несущего некоторую полезную 
информацию. 


Литература 


1. Хинтон Д.Е. Как обучаются нейронные сети / Д.Е. Хинтон // В мире науки. — 1992. — № 11-12. — 
С. 103-110. 

2. Тимченко Л.И. Многоэтапная параллельно-иерархическая сеть как модель нейроподобной схемы 
вычислений / Л.И. Тимченко // Кибернетика и системный анализ. — 2000. — № 2. — С. 114-134. 

3. Тимченко Л. Конвергентн! та дивергентн! процеси в реальних та штучних нейронних мережах / 
Л.Т. Тимченко // Вйсник ВТЦ. - 1997. — № 1.- С. 5-10. 

4. Свечников С.В. Квазиимпульсно-потенциальные оптоэлектронные элементы и устройства логико- 
временного типа / Свечников С.В., Кожемяко В.П., Тимченко Л.И. — К. : Наукова думка, 1987. — 256 с. 

5. ТиасВепКо Геот4 1. Мефо4 Юг Тгаште оЁ а РагаПе!-Нлегагсса! МебмогК, Вазе4 оп Роршаноп 
Софт» Юг Ргосеззте оЁ Ежепаеа Газег Раз Ппазез / Геоша Т. ТиасБепко, Уипу Е. Кащает, 
Уоодутуг Р. Ко7Нетуако [и др.] // РгосееЧт$ оЁ ЗРМЕ. — 2002. — Уч]. 4790. —Р. 465-479. 


«Штучний 1нтелект» 32009 615 


лот Тимченко Л.И., Мельников В.В., Кокряцкая Н.И., Кутаев Ю.Ф., Ивасюк И.Д. 


6. Сегментащя натшвтонових зображень за ознакою зв’язаност! : монографля / [Железняк А.Л., 
Тимченко Л.1., Стасюк О.Т. та 1н.]. —К. : ДЕТУТ, 2008. - 144 с. 

7. Вазоу М.О. Газег Сопёго| оЁ Меаг Еаг Зрасе ап4 РоззШ@ез ог Ветоуа[ о# Зрасе Оеылз Нот Огбй 
ул Ехр]озуе Рвою-ПО155ос1аНоп Газег$ у Рвазе Сопазайоп / М.С. Вазоу, Е.М. ХетзКоу, У.Е. 
Кищаеу её. а|. // Ргос. ОСТГ/НРГ, 98. ЭРТЕ Зутрозиит, (3{-РеегзБаго, Кизза, 1998). — $1-РаегзБиго, 
1998. 

8. Тимченко Л.И. Система координатной привязки для нестационарных сигналов / Л.И. Тимченко, 
Ю.Ф. Кугаев, А.А. Герций [и др.]| // Известия РАН. — 2001. — № 6. - С. 886-890. — (Серия 
«Физическая»). 

9. ТиасвепКо Г.1. Мефо4 Юг Ргосеззте оЁ Ежепаеа Газег Раз Ппазез / Г.Т. ТипсБепКо, У.Р. 
Ко7ретуаКо, У.Е. Клиаеу [и др.| // Адуапсез ш Еесёлса| ап Сотрщег Епетеегте — «ЗеЁп се 
Маге» / ОшуегзИу оЁ бисеауа Коташа. — 2003. — Уо/. 3(10), № 2(20). -Р. 66-78. 

10. Вступ в алгоритмачну теор!ю 1ерархй 1 паралелзму нейропод1бних обчислювальних середовищ та Й 
застосування до перетворення зображень. Основи теорй шрам1дально-стьового перетворення 
зображень / [Кожем’яко В.П., Тимченко Л.Т., Кутаев Ю.Ф., Гвасюк [.Д.]. — К. : УМК ВО, 1994. — 272 с. 

11. Авоагууа1 7.К. ОБесе гесостоп ап4 регРюгтапсе Боипа$ / Г.К. Аввагу’а1, ЗызЬи ЗВай // Ргос. Птазе 
Апа|уз15 апа Ргосеззте, (ЕЙогепсе, Па]у, Зеретфег 17 — 19, 1997). —Р. 343-360. 

12. Хьюбел Д. Глаз, мозг, зрение / Хьюбел Д. — М. : Мир, 1990. - 239 с. 

13.Рабинович З.Л. Представление и обработка знаний во взаимодействии сенсорной и языковой 
нейросистем человека / 3.Л. Рабинович, Г.С. Воронков // Кибернетика и системный анализ. — 1998. — 
№ 2. С. 3-11. 

14. [Электронный ресурс] Режим доступа : ВИр://\улу\.Коггезроп4еп. пе. 

15. Ковзель М.О. Паралельно-1ерархчне перетворення 1 О-обробка 1нформацй для систем реального 
часу : монография / [Ковзель М.О., Тимченко Л.1., Кутаев Ю.Ф. та 1н.]. - Кив : КУЕТТ, 2006. - 492 с. 

16. Тимченко Л.Г. Компактний опис моделей зображень для класиф1каци образав / Л.[. Тимченко, 
С.В. Чепорнюк, Ю.Ф. Кутаев, О.А. Герщй // Вйсник ВЦ. - 1998. — № 2. - С. 72-83. 

17. Кожем’яко В.П. Паралельно-1ерарх1чн! мереж! як структурно-функцональний базис для побудови 
спещаллзованих моделей образного комп’ютера : монографля / Кожем’яко В.П., Тимченко Л.1., 
Яровий А.А. — Ванниця : Унверсум — Вйнниця, 2005. - 161 с. 

18. ТнасБепКо Г.Г. А Вгат-ГАКе Арргоась №0 Мишзазе НегагсШа| Ппазе. Зргтеег-УеПас Ргосеззше / 
Г.[. ТиасБепко, У.Е. Каеу, 5.У. СБерогпучК [и др.] / Ргос. Птае Апа|у$15 ап4 Ргосеззте, (Е1огепсе, 
Ца[у, Зеретбег 17 — 19, 1997). -Р. 246-253. 

19. НеБь О.О. ТБе Оггатхайоп оЁ Веваутог / НеБЬ Р.О. - \Пьу, 1949. 


ЛА. Тимченко, В.В. Мельшков, Н.. Кокряцька, Ю.Ф. Кутасв, Г.Д. Твасюк 

Метод паралельно-1ераржчного розшзнавання зображень 

У стати! розглянуто метод паралельно-ерарх!чного перетворення на основ! популящиного кодування та 
його застосування для задач розшзнавання образв. Паралельно-1ерарх!чне перетворення описуеться як 
системна модель для розшзнавання образв. Наведен! теоретичн! в1домостт, експериментальн! дослдження 1 
програмна реалзащя. 
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